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Surface  discharge  in  air  has  been  studied  using  chemistry  of  formation  of  different  ion  and  neutral  species  of 
nitrogen  and  oxygen.  Continuity  equations  governing  densities  of  electrons,  ions  and  neutral  species  of  nitrogen 
and  oxygen  are  solved  with  Poisson’s  equation  to  obtain  spatial  and  temporal  profiles  of  densities,  and  voltage. 
Results  show  that  density  profiles  of  charged  species  are  nearly  opposite  to  each  other  about  the  electrode  set  at 
positive  and  negative  peaks  of  the  cycle.  Density  profiles  of  nitrogen  atom  N  and  oxygen  atom  O  are  similar  to 
each  other  at  the  same  temporal  location  in  a  cycle  because  these  do  not  respond  to  the  electric  field.  Also  V2  and 
°2  respond  to  electric  field  in  a  similar  manner,  hence,  their  density  profiles  are  similar  to  each  other  at  same 
time  instance.  There  is  difference  in  the  level  of  ionization  of  different  species  due  to  difference  in  rate 
coefficients.  The  motion  of  electrons  and  ions  in  response  to  electric  field  results  in  charge  separation.  The  value 
of  separated  space  charge  shows  footprints  of  electron  deposition  downstream  of  the  overlap  region  of  the 
electrodes  resulting  in  formation  of  virtual  negative  electrode.  The  electrodynamic  body  force  per  unit  volume 
has  been  predicted  about  the  dielectric  surface  for  the  twentieth  cycle.  A  time  averaged  force,  predominantly 
downstream  with  a  transverse  component  towards  the  wall,  acts  on  the  plasma  in  forward  direction.  This  results 
in  a  moving  wave  of  plasma  over  dielectric  surface  in  the  positive  x-direction  which  can  induce  a  control 
mechanism  in  the  neighboring  air  flow. 


Nomenclature 


D 

Diffusion  coefficient,  cm2/s 

O 

Oxygen  atom 

E 

Electric  field,  V/cm 

02+ 

Positive  oxygen  ion 

e 

Electron  charge,  statcoulomb 

O' 

Negative  oxygen  ion 

e 

Permittivity,  farad/m 

(O 

Applied  frequency,  radians 

F 

Electric  force  density,  dynes/cm3,  pN/cm3 

P 

Pressure,  Pa 

<l> 

Potential,  V 

q 

Charge  separation,  statcoulomb 

K 

Rate  coefficients,  cm3  s"1 

T 

Time  coordinate,  s 

kB 

Boltzmann  constant 

t 

Temperature,  eV 

L 

Electrode  length,  cm 

V 

Velocity,  m/s 

Ld 

Electron  Debye  length,  cm 

X 

Spatial  coordinate,  cm 

M 

mass,  kg 

P 

Mobility,  cm2V"1s"1 

Subscript 

o  * 

N 

Number  density,  cm"3 

o. 

a 

Species 

n2 

Nitrogen  molecule 

B 

Boltzmann 

N 

Nitrogen  atom 

E 

Electron,  finite  element 

N+2 

Positive  nitrogen  ion 

I 

Ion 

o2 

Oxygen  molecule 

N 

Neutral 
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I.  Introduction 

PLASMA  actuators  are  useful  flow  control  devices  for  various  aerodynamic  applications.  They  have  been 
demonstrated  for  active  mitigation  of  performance  deficiencies  under  off-design  conditions.1"4  These  actuators  are 
distinctive  for  their  absence  of  moving  parts,  rapid  on-off  deployment  and  attractive  self-limiting  characteristics.  The 
interest  of  this  work  is  in  dielectric  barrier  discharge  (DBD)  driven  plasma  actuators  operating  at  atmospheric  pressures. 
DBD  is  widely  used  in  applications  for  ozone  production,  surface  etching,  decontamination,  plasma  display  panels, 
etc.2’3  Typically  in  an  atmospheric  DBD,  a  few  watts  of  power  is  supplied  by  a  1-50  kilovolts  rms  with  50  Hz  to  50  kHz 
frequencies.  The  discharge  is  characterized  as  transient  with  a  minute  time  scale,  several  orders  of  magnitude  smaller 
than  that  for  the  bulk  gas  flow,  thus  allowing  fast  response  of  the  fluid  to  the  electrical  impulse.  The  plasma  actuators 
typically  operate  in  air.  While  successful  experiments1'4  and  simplified  numerical  modeling  efforts5'8  have  prompted  a 
tremendous  interest  in  plasma  actuators,  basic  first  principles  understanding  for  these  devices  are  yet  to  be  explored  for 
air  chemistry  so  that  the  operational  limits  can  be  further  alleviated. 

The  applications  of  plasmas  in  hypersonic  flows  for  sustained  flight  in  low-pressure  at  upper  atmosphere  and 
during  atmospheric  reentry  are  known  for  physical  complexity.9  Severe  thermal,  chemical,  and  ionizational 
nonequilibrium  conditions  are  experienced  due  to  high  velocities  and  high  altitudes  during  tranatmospheric  missions. 
The  extreme  environmental  conditions  have  a  major  impact  on  the  design  and  study  of  the  vehicle  aerodynamics  and 
thermal  loading.  Chemical  kinetic  processes  in  the  plasma  play  an  important  role  in  electromagnetic  interference 
surrounding  the  craft.  Knowledge  of  chemical  kinetic  processes  is  crucial  to  mission  planners.  The  computational  study 
has  to  include  mass,  momentum,  and  energy  equations  of  the  neutral  species  and  charge  particle  species.  The  influence 
of  the  charge  particles  on  the  flow  field  needs  to  be  included  and  the  elastic  and  inelastic  collisional  energy  transfer 
mechanisms  also  needs  to  be  taken  into  account  for  the  high  degree  of  thermal,  chemical  and  ionizational 
nonequilibrium  in  the  flow  field.  The  reactions  include  momentum  transfer,  molecular  ionization,  excitation, 
dissociation,  dissociative  recombination  and  secondary  electron  emission. 10 

The  chemistry  for  ionized  air  at  atmospheric  pressure  is  even  more  complicated  due  the  faster  species  reaction 
rates.10  Specifically,  barrier  discharges  in  air-like  N2/02  mixtures  consist  of  a  multitude  of  micro  and  nanoscale 
discharges.  Theoretical  modeling  of  these  discharge  processes  requires  a  time-resolution  in  the  subnanosecond  range  to 
describe  disparate  electron  and  ion  time  scales.  Also,  small  spatial  dimensions  of  barrier  discharges  at  these  pressures 
require  a  spatial  resolution  much  finer  than  the  submillimeter  scale.  Thus,  investigation  of  the  spatio-temporal  evolution 
of  such  multiscale  barrier  discharge  is  quite  a  challenging  but  necessary  task  to  understand  the  governing  physics  of 
discharge  in  air  and  in  N2/02  mixtures.  Barrier  discharges  offer  a  good  technical  challenge  to  properly  understand  the 
flow  and  physical  models.  Oxygen  and  nitric  oxide  species  in  the  air  could  significantly  alter  the  discharge  due  to  the 
presence  of  negative  oxygen  ions  and  nitric  oxide  species.  The  negative  ions  add  momentum  in  the  opposite  direction  to 
the  positive  ions  and  oxygen  and  nitric  oxide  have  different  recombination,  dissociation,  and  excitation  behavior.  In 
order  to  explore  the  force  production  differences  between  nitrogen,  oxygen  and  nitric  oxide,  computational  study  needs 
to  be  carried  out  by  including  major  species. 

The  asymmetric  surface  dielectric  barrier  discharge  has  been  used  in  an  experimental  setup  for  airflow 
production  close  to  the  dielectric  surface.  Effect  of  several  discharge  parameters  such  as  applied  voltage  waveform, 
distance  between  electrodes,  dielectric  thickness  and  permittivity  has  been  experimentally  determined  for  atmospheric 
pressures.  The  time-averaged  flow  velocity  spatial  profiles  has  found  the  wall  jet  velocities  of  up  to  3.5  m/s  at  heights  of 
1-2  mm.4  In  another  effort11,  the  measurements  of  the  instantaneous  flow  velocity  induced  by  surface  plasma  actuators 
in  air  at  atmospheric  pressure  have  been  conducted  with  two  different  types  of  plasma  actuators  in  order  to  determine  the 
establishment  time  of  the  induced  airflow.  A  strong  correlation  was  found  between  the  transient  regime  of  the  discharge 
current  and  the  ionic  wind  velocity.  The  AC  discharge  was  found  to  induce  a  pulsed  airflow  at  a  frequency 
corresponding  to  the  high  voltage  waveform  frequency  and  the  DC  discharge  airflow  with  important  fluctuations.11  The 
performance  of  plasma  actuators  has  been  studied  with  wind  tunnels,  drag  balances,  Pitot  tubes,  smoke  flow 
visualization,  and  fluid  dynamic  modeling  programs.  The  electric  power  transfers  the  momentum  to  the  neutral  gas  flow 
through  ion-neutral  collisions.  The  effects  of  actuator  geometry,  materials,  the  RF  frequency  and  RMS  voltage  have 
been  studied  both  experimentally3  and  numerically.12 

The  first  principles  modeling  of  the  surface  DBD  has  been  reported  in  recent  literature.  A  two-dimensional 
fluid  model  of  the  surface  discharge  showed  the  force  per  unit  volume  acting  on  the  flow  due  to  the  momentum  transfer 
from  charged  particles  to  neutral  molecules.1’  The  force  in  the  DBD  was  found  to  be  of  the  same  nature  as  the  electric 
wind  in  a  corona  discharge  and  localized  in  the  cathode  sheath  region  of  the  discharge  expanding  along  the  dielectric 
surface.  The  present  authors  have  carried  out  simulation  study  of  the  dielectric  barrier  plasma  actuator5  and  its  effect  on 
the  mitigation  of  separated  flows  over  flat  plate6’12  and  NACA0015  airfoil14"16.  They  have  shown  that  electron  deposition 
downstream  of  the  overlap  region  of  the  electrodes  results  in  formation  of  virtual  negative  electrode.  The  force  due  to 
the  charge  separation  is  in  the  positive  streamwise  direction  during  most  of  the  positive  part  of  the  rf  cycle  and  in 
negative-x  direction  during  most  of  the  negative  part  of  the  cycle.  The  magnitude  of  the  former  force  is  higher  that  the 
latter.  Thus,  domain  integration  of  the  force  on  the  charge  separation  and  time  average  of  the  force  was  found  primarily 
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to  be  in  the  positive-x  and  negative-y  direction.  In  this  paper,  we  study  atmospheric  air  discharge  using  a  mixture  of 
oxygen  and  nitrogen  about  an  asymmetric  dielectric  barrier  plasma  actuator.  We  solve  the  equations  governing 
dynamics  of  electrons,  ions  and  neutral  nitrogen  and  oxygen  species  to  obtain  spatio-temporal  profiles  of  electron 
density,  ion  density,  electric  potential,  and  neutral  gas  density.  We  identify  reasons  for  difference  in  profiles  of  different 
species.  The  value  of  separated  space  charge  and  electrodynamic  body  force  per  unit  volume  has  also  been  obtained  for 
the  twentieth  rf  cycle. 

This  paper  is  organized  as  follows.  Section  II  describes  the  problem  statement  along  with  the  air  chemistry. 
Section  III  presents  the  computational  results  and  relevant  discussion.  Section  IV  highlights  the  conclusions. 

II.  Problem  Statement 

Fig.  1  shows  schematic  of  an  asymmetric  single  dielectric  barrier  plasma  actuator.  It  consists  of  two  electrodes  separated 
by  a  dielectric.  The  upper  electrode  is  exposed  to  the  air  while  the  lower  electrode  is  placed  underneath  the  dielectric. 
The  region  simulated  is  6  cm  long  and  4  cm  high.  The  lower  part  of  the  domain  consists  of  a  0.1  cm  thick  insulator  with 
dielectric  constant  sd  =  as0  while  the  upper  part  is  filled  with  air  of^y  =  1.0055£0  ,  where  s0  is  permittivity  of  the  free 

space.  We  have  taken  a  =3.5  for  Kapton  polyimide  and  negligible  rf  electrode  thickness.  The  thickness  of  the  electrodes 
is  assumed  to  be  infinitesimally  small.  The  length  L  of  both  the  electrodes  is  2.0  mm.  The  rf  electrode  extends  from 
x=2.8  cm  to  x=3  cm  at  y=0.1  cm,  the  grounded  electrode  is  from  3.2  cm  to  3.4  cm  at  y=0,  with  a  0.2  cm  gap  between 
electrodes  along  the  x-axis.  The  embedded  electrode  is  grounded  and  a  sinusoidal  voltage  (/>  =  (/>0  sin(2 7ft)  is  applied  to 
the  exposed  electrode.  The  frequency  of  excitation  is  fixed  at  /  =5  kHz  and  excitation  potential  is  (Zi0  =  1000  volt. 


air 


4>  (0 


Dielectric 


m 

asma 


To  Ground 


FIG.  1.  Schematic  of  asymmetric  single  dielectric  barrier  plasma  actuators  operating  in  air. 


The  equations  governing  chemistry  of  discharge  are  as  follows.10  The  metastable  species  along  with  Nd  and  O4  are 
neglected  due  to  their  fast  recombination  rates.  Also,  NO  is  not  included  to  reduce  the  numerical  complexity  in  the  first 
principles  model  at  this  stage. 


(1)  Nitrogen  model 

(a)  Ionization  e  +  N2  — >  N2  +  2e , 

(1) 

(b)  Dissociation  e  +  N2  —>  2 N  +  e  , 

(2) 

(c)  Dissociative  recombination  e  +  N2  — >  N2  , 

(3) 

(2)  Oxygen  model 

(a)  Ionization  e  +  02  —>  02  +2e , 

(4) 

(b)  Dissociation  e  +  02  — »  20  +  e  , 

(5) 

(c)  Dissociative  attachment  e  +  02  —>  0_  +  O  , 

(6) 

(d)  Recombination  0~  +  02  — »  O  +  02 , 

(7) 

(e)  Dissociative  recombination  O  +  O  —>  02  +  e  , 

(8) 

The  drift-diffusion  form  of  continuity  and  Poisson’s  equations  for  the  electrons,  ions  and  neutrals  are  as  follows: 

dn 

-^-  +  V.(ne\e)  =  kxnenNi  ~k3nen^  +  (k5-k2)nen^  +k9nenQ , ,  (9) 
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Sl'l 

w  +  S7.(nNxN)  =  2k2nenN  , 
at 

(10) 

dnN 

dt  +v-(«v2Vjv2)=  kxnenNl  k2ne^N2+k3nenN, 

(11) 

dnNi 

dt  +  V-(nN;y'N;)  =  kinenN2  k3nenN+ 

(12) 

dn0 

d~  +V.(«02v02)=  (k4  +  k5+k6)n02ne+k7n0-n0+ +ksnenQ+ 

(13) 

d°  +  ^-(wovo)  =  (2&5  +  k6)nen0l  +  k1nQ  n()+ 

(14) 

dn 

dt  +V.(vvo_)  =  Mc«o2  kino  noi 

(15) 

8n0+ 

gt  +  VOj+  )  =  k4nen02  M0-'V+  k*neno; 

(16) 

sV2</>  =  4 Ke(ne  +  nQ  -  «r  -  nQ+ ) , 

(17) 

with  momentum  flux  na\a  —  —  sgn  (e) najUaV (/)  —  DaVna ,  and  electric  field E  =  -V^ . 

Above,  k] . k-x  arc  various  rate  coefficients  related  to  Equations  (1)  through  (8)  obtained  from  Kossyi  et  al.ln. 

Electron  temperature  appearing  in  these  rate  coefficients  is  calculated  from  E  =  kBTe/'Vne/ne,  which  is  obtained 
assuming  initial  Boltzmann  distribution,  n  oz  exp (e0/kBT  )  .  In  Eqns.  (9)-(  1 7)  na  and  va  are  density  and  velocities  of 
species  a.  The  mobilities  na  and  diffusion  rates  Da  are  taken  from  Ellis  et  al.17.  The  bulk  density  of  the  air  is  taken  to  be 

1.3 x  1CT3 gm/cm3.  Atmospheric  ratio  of  3.6  is  taken  for  nitrogen  to  oxygen  gas  molecules.  The  self-consistent 
formulation  is  solved  using  a  Galerkin  variational  formulation  based  finite-element  method18  to  obtain  electron,  ion  and 
neutral  species  densities  of  nitrogen  and  oxygen,  and  the  electric  potential. 

Initial  and  boundary  conditions  are  as  follows.  The  all  initial  particle  concentrations,  except  those  of  the 
electrons,  nitrogen  and  oxygen  molecules,  are  taken  to  be  zero.  Initial  electron  density  is  taken  lOVcm3.  Initial  rf 
potential  is  also  zero.  All  neutral  particle  concentrations  are  taken  zero  at  the  boundaries  of  the  dielectric  and  normal 
components  of  flux  are  zero  at  outer  boundaries  of  plasma.  Normal  components  of  electric  field  are  taken  zero  on  all  the 
outer  boundaries  of  the  domain  (dielectric  as  well  as  plasma  domain).  No  electron  and  ion  current  flow  inside  the 
dielectric  is  considered.  Electric  field  normal  to  the  dielectric  surface  is  discontinuous  by  the  separated  charge. 


III.  Results  and  Discussion 

As  we  switch  on  the  rf  voltage,  discharge  in  the  air  occurs.  Different  ion  and  neutral  species  are  formed  through 
ionization,  dissociation  and  dissociative  attachment  etc.  Recombination  also  occurs.  Full  chemistry  of  different  species 
formation  and  recombination  is  described  by  Equations  (1)  to  (8).  Figure  2  (a)  shows  electron  density  ne  as  a  function  of 
normalized  time  cot  for  x/L=15  and  y/  /,H  =25.  Steady  state  is  reached  within  a  few  cycles,  thereafter;  there  is  a  small  loss 
in  the  densities,  as  described  in  Figure  2a. 

Initial  electron  density  is  taken  10  '  /cm3.  The  electrons  are  repelled  from  the  area  of  grounded  electrode  and  are 
attracted  towards  rf  electrode  during  positive  part  of  the  cycle.  Reverse  happens  during  negative  part  of  the  cycle,  hence, 
density  profiles  are  opposite  of  each  other  in  Figure  2b-c  in  the  context  of  location  of  electrodes.  Spatial  and  temporal 
density  profile  of  electrons  is  governed  by  continuity  Equation  (10)  coupled  with  other  equations.  Electron  density 
grows  to  a  level  1011  /cm3.  Electrons  respond  to  electric  field  very  fast  due  to  low  mass  and  high  mobility.  There  is 
deposition  of  electrons  on  the  dielectric  surface  over  grounded  electrode  during  negative  part  of  the  cycle.  Dielectric 
surface  behaves  as  a  negative  electrode  after  negative  peak  of  the  rf  voltage  till  the  end  of  the  cycle. 
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(a) 


0  10  20  30  uO  10  20  30 


x/L  x/L 

(b)  (c) 

FIG.  2.  Evolution  of  electron  density  ne  (a)  as  a  function  of  normalized  time  at  for  x/L=15  and  y /  AD  =25,  (b)  at 
positive  peak  of  the  20th  cycle  at  =  38.5/7" ,  and  (c)  at  negative  peak  of  the  cycle  at  =  39.5zr . 


Initial  density  of  the  nitrogen  atoms  was  taken  zero.  The  nitrogen  atoms  are  produced  through  dissociation. 
Chemistry  of  nitrogen  atom  formation  is  given  by  Equation  (2)  and  its  density  is  governed  by  continuity  Equation  (10) 
coupled  with  other  equations.  Density  of  nitrogen  atoms  grows  to  the  same  level  as  that  of  electrons,  i.e.,  10n  /cm3. 
Figure  3  shows  distribution  of  nitrogen  atom  density  at  positive  and  negative  peaks  of  the  twentieth  cycle.  Since 
nitrogen  atom  do  not  have  any  charge,  no  electric  force  works  on  these  atoms. 
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FIG.  3.  Distribution  of  nitrogen  atoms  N  (a)  at  positive  peak  of  the  cycle  at  =  38.5zr  ,  and  (b)  at  negative  peak  of 
the  cycle  cot  =  39.5 n . 
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Initial  density  of  the  nitrogen  ions  N2  was  taken  zero.  The  nitrogen  ions  are  produced  through  ionization 
of  nitrogen  molecules.  Chemistry  of  nitrogen  ions  formation  is  given  by  Equation  (1)  and  its  density  is  governed  by 
continuity  Equation  (12)  coupled  with  other  equations.  Density  of  nitrogen  ions  A^  grows  to  the  same  level  as  that  of 
electrons,  i.e.,  1011  /cm3.  Figure  4  shows  distribution  of  nitrogen  ions  density  n  ^ ,  at  positive  and  negative  peaks  of  the 

twentieth  cycle.  The  nitrogen  ions  A^  are  repelled  from  the  area  of  rf  electrode  and  are  attracted  towards  grounded 
electrode  during  positive  part  of  the  cycle.  Reverse  happens  during  negative  part  of  the  cycle,  hence,  density  profiles  are 
opposite  of  each  other  in  Figure  4.  Since  positive  ions  respond  to  electric  field  in  an  opposite  manner,  density  profiles 
are  opposite  to  that  of  density  profile  of  electrons  in  the  context  of  location  electrodes.  The  ions  are  much  heavier  than 
that  of  electrons;  mobility  of  ions  is  low,  the  ions  remain  close  to  the  wall  than  electrons. 
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FIG.  4.  Distribution  of  nitrogen  ions  Nl  density  (a)  at  positive  peak  of  the  cycle  cot  =  38.5zr  ,  and  (b)  at  negative 
peak  of  the  cycle  cot  =  39.5zr  . 


Initial  density  of  the  oxygen  atoms  O  was  taken  zero.  The  oxygen  atoms  are  produced  through  dissociation  of  oxygen 
molecules.  Chemistry  of  oxygen  atoms  formation  is  governed  by  Equation  (5)  and  its  density  is  governed  by  continuity 
Equation  (14)  coupled  with  other  equations.  Figure  5  shows  density  profile  of  oxygen  atom  O  at  positive  peak  of  the 
cycle  cot  =  38.5;r  ,  and  at  negative  peak  of  the  cycle  cot  =  39.5 n  .  Density  of  oxygen  atoms  is  higher  than  that  of  nitrogen 
atoms  due  to  difference  in  rate  coefficients  related  to  formation  of  nitrogen  and  oxygen  atoms.  Density  profiles  of  both 
nitrogen  and  oxygen  atoms  are  similar  because  both  are  neutral  atoms  and  do  not  respond  to  electric  field. 


(a)  (b) 


FIG.  5.  Distribution  of  oxygen  atom  O  density  (a)  at  positive  peak  of  the  cycle  cot  =  38.5zr  ,  and  (b)  at  negative 
peak  of  the  cycle  cot  =  39.5zr  . 
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Initial  density  of  the  oxygen  ions  O  was  taken  zero.  The  oxygen  ions  O  are  produced  through  dissociative 
attachment  of  oxygen  molecules  with  electrons  governed  by  Equation  (6).  Its  density  is  governed  by  continuity  Equation 
(15)  coupled  with  other  equations.  Density  of  oxygen  ions  O  grows  to  a  level  of  slightly  less  than  that  of  density  of 
nitrogen  ions  N2  ■  Figure  6  shows  distribution  of  density  of  oxygen  ions  O  at  positive  and  negative  peaks  of  the 

twentieth  cycle.  The  oxygen  ions  O  are  repelled  from  the  area  of  grounded  electrode  and  are  attracted  towards  rf 
electrode  during  positive  part  of  the  cycle.  Reverse  happens  during  negative  part  of  the  cycle,  hence,  density  profiles  are 
opposite  of  each  other  in  Figure  6.  The  oxygen  ions  O  are  much  heavier  than  that  of  electrons  and  mobility  is  low; 
response  is  delayed.  The  ions  remain  close  to  the  wall  than  electrons. 
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FIG.  6.  Distribution  of  negative  oxygen  ion  O  (a)  at  positive  peak  of  the  cycle  cot  =  38.5zr  ,  and  (b)  at  negative 
peak  of  the  cycle  cot  =  39.5zr  . 


Initial  density  of  the  oxygen  ions  02  was  taken  zero.  The  oxygen  ions  ()2  are  produced  through  ionization  of 
oxygen  molecules.  Chemistry  of  oxygen  ions  02  formation  is  given  by  Equation  (4)  and  its  density  is  governed  by 
continuity  Equation  (16)  coupled  with  other  equations.  Density  of  oxygen  ions  02  grows  to  the  same  level  as  that  of 
02  i.e.  5  x  1 0 10  /cm3.  Figure  7  shows  distribution  of  oxygen  ions  density  02  at  positive  and  negative  peaks  of  the 

twentieth  cycle.  The  N2  and  02  respond  to  electric  field  in  a  similar  manner,  hence,  their  density  profiles  are  similar 
to  each  other  at  same  time  instance.  The  difference  between  density  levels  of  nitrogen  and  oxygen  is  because  of  ratio  of 
these  gases  in  atmospheric  air. 
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FIG.  7.  Distribution  of  positive  oxygen  ion  02  (a)  at  positive  peak  of  the  cycle  cot  =  38.5zr  ,  and  (b)  at  negative 
peak  of  the  cycle  cot  =  39.5zr  . 


7 

American  Institute  of  Aeronautics  and  Astronautics  Paper  No.  2006-3381 


Charged  particles  move  in  response  to  the  electric  field  as  described  above.  This  give  rise  to  charge 
separation  ijr  =  c{n()  +  n  y  —ne—  n()  ) .  Density  profile  of  charge  separation  is  shown  in  Figure  8.  It  resembles  that  of 

density  profiles  of  positive  nitrogen  and  oxygen  ions.  Level  of  charge  separation  matches  density  of  ions  and 

electrons  because  densities  of  other  charge  species  <9  and  is  small  in  comparison.  Charge  separation  gives  rise  to  a 
self-generated  electric  field  through  Poisson’s  Equation. 


FIG.  8.  Distribution  of  charge  separation  q  =  e(nQ+  +  «w+  -  ne  -  nQ-)  (a)  at  positive  peak  of  the  cycle  cot  =  38.5^  , 
and  (b)  at  negative  peak  of  the  cycle  cot  =  39.5;r  . 


Figure  9  shows  charge  separation  q  =  e(/?0+  +  nN+  ~ne~  nQ)  along  the  surface  of  dielectric  at  y/AD=2 0  and 

y/AD=  40.  Charge  separation  is  highest  in  the  area  of  electrodes  where  discharge  activity  occurs.  Charge  separation  is 
highest  at  voltage  peaks  of  the  cycle  because  of  high  electric  field.  Low  electric  field  is  associated  with  cot  =  39/T  and 
40/r for  yl XD=  40,  hence,  charge  separation  is  low  at  these  points.  We  can  also  see  that  peak  of  charge  separation 
moves  towards  higher  value  of  yt  AD  at  positive  peak  of  the  cycle.  A  small  negative  charge  deposition  is  also  noticed. 
These  solutions  are  basic  and  needs  to  incorporate  the  effect  of  nitrous  oxide  and  metastable  species. 


(a)  (b) 

Figure  9.  Charge  separation  q  as  a  function  of  normalized  distance  x/L  for  cot  =  38.5zr ,  39;r ,  39.5;r  and  40;r  at 
(a)  y/AD=20  and  (b)  y/AD=  40. 


Electric  field  E  exerts  a  force  F  =  e(nN+  +  nQ+  -  ne  -  n()  )E  on  the  charge  separation.  Figure  10  shows  vectors 

of  force  distribution  in  the  domain  at  positive  and  negative  peaks  of  the  twentieth  cycle.  The  force  is  directed  towards 
positive-x  direction  at  positive  peak  of  the  cycle  and  towards  negative-x  direction  at  negative  peak  of  the  cycle.  It  is  also 
directed  towards  wall  at  both  the  points  of  the  cycle.  The  force  follows  direction  of  the  electric  field. 
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FIG.  10.  Vectors  of  force  in  the  domain  (a)  at  positive  peak  of  the  cycle  cot  =  38.5;r  ,  and  (b)  at  negative  peak  of 
the  cycle  cot  =  39.5;r  . 

Figure  1 1  shows  time  average  of  the  force  during  twentieth  cycle.  The  time  average  of  the  force  is  mostly 
directed  towards  negative  y-direction  with  orientation  towards  positive  x  direction.  The  force  lines  are  concentrated  over 
the  exposed  and  grounded  electrodes.  The  electrons  respond  to  the  electric  field  very  quickly  due  to  high  mobility.  The 
electrons  get  deposited  at  the  surface  of  the  dielectric  above  grounded  electrode  during  negative  part  of  the  cycle  and  the 
dielectric  surface  becomes  negatively  charged.  It  behaves  as  virtual  negative  electrode  after  negative  peak  of  the  cycle 
till  the  end  of  cycle.  The  peak  magnitude  of  the  force  is  located  downstream  of  the  rf  electrode. 


FIG.  11.  Vectors  of  time  average  of  force  distribution  in  the  domain  during  a  full  cycle,  x  and  y  are  in  cm. 

The  x-component  of  the  force  Fx  =  qEx  is  shown  in  Figure  12  along  the  surface  of  dielectric  at  y/AD=20  and 
y/AD=40.  Magnitude  of  force  is  large  at  the  peaks  of  the  cycle.  Most  of  the  force  is  concentrated  in  the  region  of 
electrodes.  The  force  is  positive  at  positive  peak  of  the  cycle  and  negative  at  negative  peak.  The  y-component  of  the 
force  is  mostly  negative  as  shown  in  Figure  13.  We  can  see  again  that  force  is  highest  at  the  peaks  of  the  cycle  and 
concentrated  in  the  region  of  electrodes.  At  positive  peak  of  cycle,  peak  of  the  force  is  at  higher  value  of  x/L.  These 
results  can  be  explained  by  considering  direction  and  magnitude  of  the  electric  field.  Direction  of  electric  field  is  from  rf 
electrode  to  grounded  electrode  at  positive  peak  of  the  cycle  and  vice  versa.  Magnitude  of  electric  field  is  highest  at  the 
peaks  of  cycle.  Electric  field  decreases  with  distance  from  electrodes,  hence,  force  also  decreases. 
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Figure  12.  x-component  of  the  force  Fx  as  a  function  of  normalized  distance  x/L  for  cot  =  38.5;r  ,  39n ,  39.5;r  and 
40 n  at  (a)  y!) 1^=20  and  (b)  y/2D  =40. 


Figure  13.  y-component  of  the  force  Fv  as  a  function  of  normalized  distance  x/L  for  cot  =  38.5/r ,  39/r ,  39.5/r  and 
40;r  at  (a)  y!) lD=20  and  (b)  y/20  =40. 


IV.  Conclusions 

Air  discharge  has  been  studied  using  an  asymmetric  dielectric  barrier  plasma  actuator.  A  self-consistent  plasma 
actuator  model  is  employed  to  study  surface  discharge  of  air.  Chemistry  of  formation  of  different  ion  and  neutral  species 
of  nitrogen  and  oxygen  has  been  taken  into  account.  Species  with  very  high  recombination  rates  and  nitrous  oxide  have 
been  neglected  for  simplicity.  Continuity  equations  governing  densities  of  electrons,  ions  and  neutral  species  of  nitrogen 
and  oxygen  are  solved  with  Poisson’s  equation  using  a  two-dimensional  finite  element  based  formulation  of  plasma  to 
obtain  spatial  and  temporal  profiles  of  densities,  and  voltage.  Poisson’s  equation  governs  electric  field  generated  due  to 
charge  separation.  Computed  results  show  that  for  this  simplistic  model  oppositely  charged  species  move  in  opposite 
directions  due  to  applied  driver  rf  potential.  Response  of  electrons  to  electric  field  is  very  fast  as  compared  to  ions  due  to 
small  mass  and  high  mobility  of  electrons.  Density  profiles  of  different  species  are  nearly  opposite  to  each  other  at 
positive  and  negative  peaks  of  the  cycle.  Density  profiles  of  nitrogen  atom  N  and  oxygen  atom  O  are  similar  because 
these  do  not  move  due  to  electric  field.  Density  profiles  of  /Vj  and  02f  are  similar  to  each  other  at  any  specific  time 
instance  because  they  respond  to  electric  field  in  a  similar  manner.  There  is,  however,  difference  in  level  of  ionization  of 
different  species  due  to  difference  in  rate  coefficients.  The  motion  of  electrons  and  ions  results  in  charge  separation. 
Electron  deposition  downstream  of  the  overlap  region  of  the  electrodes  during  negative  part  of  the  cycle  results  in 
formation  of  virtual  negative  electrode.  The  value  of  charge  separation  q  =  e( nQ  +  n  v  -  ne  -  nf)  ) ,  and  the 

electrodynamic  force  per  volume  F  =e(nN+  +  n()  -  ne  -  n()  )E  have  been  obtained  near  the  dielectric  surface  for 

twentieth  cycle.  Time  average  of  the  x-force  is  positive  and  y-force  is  negative  over  the  domain,  therefore,  there  is  an 
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average  net  force  on  the  plasma  in  positive  x  and  negative-y  direction.  This  will  result  in  a  moving  wave  of  plasma  over 
the  dielectric  surface  in  positive  x-direction  which  can  find  application  in  flow  control.  A  detail  air  chemistiy  model 
with  Poisson-Navier-Stokes  equations  is  also  under  development  for  exploration  of  the  effect  of  a  plasma  actuator 
operating  in  air  using  realistic  geometry  and  electrode  configurations.  In  this  model,  an  in-depth  study  involving  more 
species  like  NO  will  be  incorporated  to  investigate  the  effect  of  surface  air  discharge  and  its  coupling  with  the  neutral 
gas  flow. 
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